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THE USE OF UNEAR FEATURE DETECTION TO INVESTIGATE 
THEMATIC MAPPER DATA PERFORMANCE AND PROCESSING 


INTRODUCTION 


The objective of the w<»’k described is to investigate the geometric and radiometric 
characteristics of Thematic Ma(^>er data throi^ analysis of linear features in the 
data. The particular aspects considered are: 

1. Thematic Mappo* ground IFOV 

2. Radiometric contrast between linear featires and background 

3. PrecisicH) of system geometric correction 

4. Band-to^nd registration 

5. Potential utility of TM data for linear feature detection especially as 
compared to MSS data. 

MODEL DESCRIPTION 

In the context of this investigatim linear features are deflned as two close, parallel 
and opposite e<%es. Examples in remotely sensed data are features sudi as roads, 
rivers and bridges. 

These features may be detected uring a local operator which considers the 
arrangement of grey levels within a S by 5 pixel array. Details are given in Gurney 
(Ref. 1), Md the m^el is described teiefly below. Fourteen possible orientations of 
the linear feature in the local area are considered, these are illustrated in Fig. 1. 
For each orientation the sum of values in column B b eom(wred to that in columns A 
and C. For dark linear features the sum in B must be less than both the sums in A 
and C by a given tturMhold. In addition B2 must be lower in value than both A2 and 
C2. For bright linear features the sums in A and C must be less than those in B. If 
these conditions are satisfled for any one of the 14 orientatims then pixel B2 is 
considered as linear. 
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Fig. 1 Orientation of a Linear 
Feature Three Pixels Long 
Passing Through the Centre of 
an Array 


Fig. 2 Variations in Threshold with 
Linear Feative Alignment 
with the Pixel Boundaries 


Results d^>end on two factors: the contrast between line material and background 
material, and the width of the feature. Any given feature may be assumed to be of 
constant vddth and contrast over the area of study. Therefore by varying the 
detector threshold it is possible to vary the accuracy with whidi the feature Is 
detected. At low thresholds all of a feature may be detected, as the tiveshold b 
increased less and less of the feature will be detected. 

An additional factor which must be considered b the alignment of the linear feative 
with the pixel boundaries. A simplified schematic represenUtion b illustrated in 
Fig. 2. If a feature occupies only a single pixel then a higher threshold may be used 
to detect it than if it b equally divided between two pixeb. Thus there b a wide 
variability between Uwesholds required to detect the full range of conditions. 

If it b assumed that the present material in a pixel b linearly proportional to its 
contribution to the reflectance, and if it b assumed that every alignment of the 
feature with the pixel boundaries b equally prc^ble then it b posable to derive a 
set of expressions relating detector threshold to detection accuracy, and faattre 
width and contrast. Complete derivations are given in Oumey (Ref. 2). 
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For Unas up to 1 pixel wide: 
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8bW-(SbW-bW/2)A 


where 


T B detector ^eshold 
A B % eecuracy / 100 

b ~ c<ffitra8t “ ( Imean backgrcnind • mean Une I )/100 
W « Une width as % pixel width 


Similarly for Unes 1 > 2 pixels in width: 
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Similar e:q>ressions may be derived for features ip to four pixels in width. Beyond 
this value detection using this model will not be po^ble. 


Contrast measurements may be made by selecting training areas corre^nding to 
known Unear and background materials and using the mean grey level values in these 
areas. Since the contrast measurement is based iq>on mean values it is clear that a 
high variance in a training area will lead to results which e<mform less closely to the 
theoretical expectations of the modeL 


Using this model it is possible to select good thresholds for features of knovm width 
and contrast. Also, by ai^lyii^ the procedure at a givoi theshold and calculating 
observed detection accuracy then it is possible to estimate feature widths, assuming 
that the pixel size is known. If aU parameters are known then it is possible to <^tain 
an estimate of pixel size. AU of these appUcations have been tested using d^lital 
MSS and RBV data with good results. 

SCENE CONTRAST AND VARIABILITY 


Using MSS data it is generally the ease that only one, or maybe two, bands provide 
sufficient contrast to allow Unear feature detection. With TM data there are more 
bands providing good contrast. This allows the possibiUty of independent verification 
firom band to band of feature width estimations. 


Figs. 3 > 10 illustrate the results of Unear feature detection using part of a scene of 
Iowa for August 1982. Each of Figs. 3-8 corresponds to detection of twight Unear 
feature in bands 1, 2, 3, 4, 9 and 7 u^ng the same toreshold (9) in each case. 
Variations between results therefore correspond to variations in contrast and scene 
variabiUty. Band 6 data, corresponding to the thermal infra*red channel, were not 
used since no features were detected. 


Results indicate that bands 1 and 3 give the best performance, with band 2 providing 
similar results. Bands 5 and 7 also perform well, but there is evidence of increased 
^oise* i.e., pixels detected as Unear which do not eorreq>ond to any recognizable 
ground feature. 
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Fig. 9 Results of Linear Feature Fig. 10 Results of Linear Feature 
Detectiem for Iowa Test Area, Detection for Iowa Test Areay 

Band 4, Dark Lines Band 5, Dark Lines 


Figs. 9 and 10 show the results of dark feature detection using bands 4 and 5. Other 
bands were unabie to detect dark features successfully. Bai^ 5 P«*fwms t>«t in 
isolatirc the rivers without picking vp excessive amounts of noisj^ Band 4 alM tww 
wi interesting effect of evenly spac^ ’dots’ along field edges. The reasons for this 
effect are not known at present and will receive further study. 


ESTIMATION OF LINEAR FEATURE WIDTHS 


Estimation of feature widths using TM data was carried out in order to establish a 
lower limit on the width of detectable feature, and to confirm that the data are 
suitable for application of the model described previomly. 


Two test areas were selected which had different ground characterises. The flrrt 
oorresponds to an area in North Carolina with a very homogeneous, fiat background 
traversed by narrow gravel surfaced roads less than a pixel in width. The second 
corresponds to part of the Washli^ton, D.C. beltway with a width between one and 
two pixels and traversing a much more variable background. 
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For Mch test aroa trainii^ areas were aelected to eatablish line: background 

contrast. The line detection procedure was then applied at at least three different 
thresholds to each band of data for which cmtrast was sufficient to allow some 
degree of feature detection. For each of these the resulting detection accuracy was 
then calculated. Results were substituted in equations 1 and 2 in order to obtain an 
estimate of the feature width as a percent of the pixel width. 


Tables 1 and 2 and Fig. 11 show results for the North Carolina area. 


Table 1. 

UNE: BACKGROUND CONTRASTS FOR NORTH CAROLINA TEST AREA 


Band 

Number 

Backeround 


Line 

Contrast 

Deviation 

1 

X 

64.5 

a 

1.5 

X 

107.4 

o 

14.9 

.43 

38.4 

2 

24.6 

0.9 

54.6 

9.9 

.30 

36.3 

3 

21.7 

2.6 

69.3 

15.6 

.48 

38.3 

4 

62.7 

4.0 

81.8 

10.4 

.19 

75.4 

5 

41.8 

6.8 

139.6 

26.4 

.98 

33.9 

6 

— 

— — 

— — 

II i 

— 

— 

7 

12.4 

2.8 

79.0 

21.0 

.67 

35.8 


Contrast - (Imean background - mean line|) /lOO. 
Deviation ^ (abackground * linea)/contrast. 


Table 1 gives the contrast data and an additional parameter termed the deviation. 
The deviation is dependent on the amount of variatim in the line and background 
data, and it essentially indicates the degree to which results are likely to accord with 
model predictions. High variances coupled with low contrast leads to considerable 
deviation from expected results. From previojs experience using MK data devl^^ion 
values of 40 or less indicate very close cmespondence between observed and 
predicted results. 
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Table 2. 

RESULTS FOR NORTH CAROLINA TEST AREA. 



Detector 

Observed Accuracy at 

Estimated Road Wi 


Threshold 

95% Confidence Level (%) 

as % Pixel Size 

BANDl 

9 

81.7 - 87.5 

21.8 - 25.8 


15 

55.5 - 63.3 

21.6 > 24.6 


18 

35.4 - 43.2 

19.8 - 21.8 


21 

25.4 - 32.6 

20.4 - 22.1 

BAND 2 

9 

37.9 - 45.7 

14.6 - 16.2 


12 

32.5 - 40.1 

18.3 - 20.0 


15 

16.6 - 23.0 

19.3 - 20.6 


18 

7.1 - 11.7 

21.2 - 22.2 

BANDS 

12 

75.0 - 81.6 

22.2 - 26.0 


15 

51.6 - 59.4 

18.3 - 20.6 


18 

36.8 - 44.6 

18.0 - 19.9 


21 

26.1 - 33.3 

18.6 - 20.2 

BAND 4 

9 

50.8 - 58.6 

27.3 - 30.9 


12 

34.4 - 42.0 

29.5 > 32.4 


18 

13.0 - 18.8 

35.7 - 37.4 

BANDS 

21 

69.3 • 76.3 

16.9 > 19.6 


24 

63.0 - 70.4 

17.2 - 19.7 


27 

57.4 - 65.2 

17.6 - 20.1 


36 

45.1 - 53.1 

19.6 - 22.0 

BAND? 

18 

61.1 - 68.7 

18.2 - 20.9 


21 

48.4 > 56.4 

17.5 - 19.7 


24 

42.2 - 50.2 

18.4 - 20.5 


27 

34.8 • 42.6 

18.9 - 20.8 



Mean road width (%) 

JS. 

Road width for 
28.5 pixel 

Excluding Band 4 

20.0 

2.0 

5.7 metres 

Including Band 4 

21.5 

4.7 

6.1 metres 


Table 2 lists the observed accuracies and calculated range of road widths for all 
band^ and thresholds used. Resulis show a close correspondence between estimated 
widths at different bands and thresholds, with the exception of results for Band 4. 
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OF POOR QUALITY 


Btnd 4 glvet a largar main width and also shows mort variability in tha astimatas for 
diffarent thrasholds. This is to ba axpaotad from tha highar daviatlon valua obsarvad 
for this band. Tha astimatad road width from all thasa data is about 6 matras 
aaauming a S8.5 by S8.S matra pixaL 



Fig. 11 Predicted vs. Observed Results for N. Carolina Test Area 


Pig. 11 shows a plot of predicted results assumii^ a 6 metre road width, against 
actual observed results. There a very close correspondence for all but the Band 4 
data. The actual width of the roads in this area is not known at the time of writing, 
but is believed to be close to the estimate giv«i hare. 

Table 3. 

LINE: BACKGROUND CONTRASTS FOR WASHINGTON, D.C. TEST AREA 


Band 

Number 

Backcround 


Line 

Contrast 

Deviation 

1 

X 

56.6 

c 

1.6 

X 

73.2 

o 

8.5 

.17 

59.4 

2 

21.6 

0.8 

62.1 

5.0 

.11 

52.7 

3 

22.8 

1.06 

66.6 

7.3 

.14 

66.1 

4 

66.4 

2.96 

43.1 

2.4 

.07 

76.6 

5 

68.6 

5.7 

60.7 

10.1 

.22 

71.8 

6 

llS.l 

0.8 

115.7 

0.5 

.01 

130.0 

7 

13.1 

2J 

28.2 

6.4 

.15 

71.3 


Contrast ■ ( | Mean background - Mean line | )/100 
Deviation « (obackground c Une)/eontrast. 


Tftblet 3 and 4 and Figi. and 13 show rasults for the Washii^on, D.C* araa. Band 
6 and Band 4 ware included in the contrast measurementt hut were not used 
•Uhsequently because of the poor contrast levels oteerved. In general, contrasts are 
considerably lower than for the North Carolina area reflectir^ both the difference in 
materials and the lower sun angle of this November scene. Deviation values are also 
correspondingly higher. 


Table 4. 

RESULTS FOR WASHINGTON, D.C. TEST AREA 


Detector Observed Accuracy % Estimated road width 

Threshold at 95% Confidence level m % pixel ria e 


BAND 1 

15 

83.6 - 94.6 

172.5 - 176.9 


18 

69.9 - 84.7 

170.1 - 176,2 


21 

47.9 - 65.3 

171.2 - 179.5 


24 

45.9 - 63.5 

164 J - 172.6 


BAND 2 

9 

65.5 > 81.2 

180.3 - 187.8 


15 

31.9 - 49.3 

173.1 - 182.1 


18 

19.5 - 35.1 

167.3 - 176.1 


BAND 3 

9 

87.4 - 97.0 

179.8 • 183.9 


15 

68.3 - 83.3 

170.1 - 176.4 


13 

58.9 • 75.5 

165.1 - 171.9 


BAND 5 

9 

58.0 > 74.8 

198.8 - 209.3 


18 

32.7 - 50.1 

196.9 - 209.6 


27 

17.3 - 32.7 

189.1 - 200.8 

BAND 7 

9 

75.2 - 88.8 

184.7 - 191.3 


15 

31.9 - 49.3 

189.3 - 201.1 


18 

9.1 - 22.1 

198.3 - 210.0 



Mean road width (%) 

£ 

Road width for 
28.5 m pixel 

Bands 1-3 

174.96 

5.14 

49.9 metres 

Bands 5, 7 

198.30 

6.63 

56.5 me^M 
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Estimated road widths are given in Table 4. These ^w a close correspondence 
between result'^ in bands 1, 2, and 3 giving a width estimate of 175% of a pixel which 
is close to SO metres for a 28.5 m pixel. The true width of the beltway is known to 
be 48-50 metres in this area so that this result corresponds very closely to the true 
width. However results for bands 5 and 7 give a width estimate close to 56 metres 
which is an over-Mtimate of the width. 




Figs. 12 and 13 illustrate {M'edicted against observed results assuming a road with of 
50 metres. The correspondence is less close than for the North Carolina data, and is 
very poor indeed for bands 5 and 7. The deviation values for these bands are both 
very high. However, it is also possible that the Band 7 result is due to a mis- 
estimation of the contrast since the slope of the observed line could intercept tho y 
axis at the same location as all other bands^ % required by the modeL 

It is concluded that the model developed here may be used to estimate the width of 
unknown features in TM data. In optimum conditions with high contrast and 
homogeneous materials it appears possible to measure widths as little as one-quarter 
of a pixeL This has considerable potential for mapping purposes. Contrast levels 
observed compare favorably with those for MSS data, although variances and hence 
deviation levels tend to be greater. However, the increased dimensionality of the 
TM data allows independent verification of results from band to band so that the 
results may be considered to be more reliable. Pixel sise Is confirmed by these 
experiments to be 28.5 metrM. 

BANDt BAND REGISTRATION 

Using linear feature detection results it is also possib'e to investigate band to band 
registrati<m. Registratim differences will be apparent as shifts in the locations of 
detected features from band to band. 


OF POOR QUALITY 



Fig. 14 1-95 Interchange, Bands 1 and 

3 Superimposed 


Fig. 15 1-95 Interchange, Bands 1 and 

7 Superimposed 


Figs. 14 and 15 illustrate the superpositicm of results of linear feature detection from 
part of the Washii^^on, D.C. scene correspondir^ to the 1-95: Beltway interchange. 
Fig. 14 shows bands 1 and 3 together; band 1 results are in black and band 3 in gray. 
There is clearly a close cwrespondence between the two. Fig. 15 shows bands 1 and 
7 together. In this case there is a clear mis-registration with the band 7 data being 
shifted both to the right and down with respect to the band 1 data. This result was 
also observed for band 5 data against band 1. Bands 5 and 7 showed a slight 
difference that was difficult to assess because of the poor detection rate in band 5. 

This method is less suitable for measurement of misregistration than it b for 
graphically illustrating that it exists. However, it does appear that there b a shift 
between primary and cold focal planes of less than 1 pixel in both directions that b 
apparently greater In the x direction than in the y directim. 

CONCLUSIONS 

It has been shown that TM data may be used to estimate TM pixel rize and to 
illustrate band: band mis-registration. Further, the geometry and radiometry of the 
data are sufficiently i»'ecbe to allow accurate estimation of the widths of linear 
features. In optimum conditions features one quarter of a pixel in width may be 
accurately measured. These results have considerable potential for applications for 
hydrological and topographic mapping. 
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